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We search for decays of a -B meson into a neutral D meson and a charged kaon, with the D 
meson decaying into a charged kaon, a charged pion, and a neutral pion. This final state can be 
reached through the 6 — > c transition B~ — > D°K~ followed by the doubly-Cabibbo-suppressed 
D° -> K+TT~Tr°, or the 6 m transition S" -> D°K- followed by the Cabibbo-favored D° 

K'^-K~-K^ . The interference of these two amplitudes is sensitive to the angle 7 of the unitarity 
triangle. We present results based on 226x10® e+e~ — > T(AS) — > BB events collected with the 
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BABAR detector at SLAC. We find no significant evidence for these decays and we set a limit 



Ral 



r([ff+jr~7r"]pK~) + r([X~7r+7r°]£)i<'+) 



r([if + ^r-7^0]Di^+) + ^{[X-7^ + 7rO]oif- 
with a Bayesian approach into rs = \A{B' 
level. 



< 0.039 at 95% confidence level, which we translate 
D°K-)/A{B- D°K-)\ < 0.19 at 95% confidence 



PACS numbers: 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 



measure ratios of decay rates: 



Following the discovery of CP violation in B meson 
decays and the measurement of the angle f3 of the unitar- 
ity triangle 7| associated with the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix, the focus has 
turned toward the measurements of the other angles a 
and 7. Following Ref. 0, several methods have been 
proposed to measure the relative weak phase between 
the D^K^ amplitude, proportional to the CKM 

matrix element Vcb (Fig. [T|), and the B~ D^K~ am- 
plitude, proportional to Vub- This weak phase, which by 
definition is 7 = arg(— can be measured 
from the interference that occurs when the D° and the 
decay to common final states. 




B- 



\1. 




FIG. 1: Feynman diagrams for the CKM-favored B — 
D°K~ and the CKM- and color-suppressed B~ D°K' 
decays. 



As an extension of the method proposed in Ref. Q , we 
search for B- [K+n-n°]DK- [|, where the CKM- 
favored B^ D^K~ decay, followed by the doubly 
Cabibbo-suppressed 13° — > K^tt^tt^ decay, interferes 
with the CKM-suppressed B^ D'^K^ decay, followed 



by the Cabibbo-favored D" ^ 



decay. 



In order to reduce the systematic uncertainties, we 
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Rads = 



T{[K+tt-tt°]dK-)+T{[K-tt+tt^]dK+) 



A 



ADS 



= rl+rl) + 2rBrDC cos J, (1) 



T{[K+7r-TrO]DK-) + T{[K- 
= 2rBrDSsm'y/RADS, 



where rs 



A{B-^D°K~) 
A{B-^D»K-) 



2 _ BiD^^K+Tv-Tv") 
'^D — B{D"^K--K + Tr») ■ 



(2) 



The 



C and S parameters are defined as: 



_ / AD{m)AD{rn) cos((5(to) - 6{rn) + SB)drn 
O — = , (6) 



, (4) 



'/ |,4£)(m)pfiTO • J |^£)(TO)|2(im 
/ AD{rn)AD{rn) sin(^(TO) — S{rn) + SB)drn 



where m indicates a point in the Dalitz plane 
['^Ktt^'^k-^o], [AD{m),S{m)] {[AD{m),6{rn)])^ the ab- 
solute value and the strong phase of the D° (-0°) decay 
amplitude, and Sb the strong phase difference between 
the two interfering B decay amplitudes. Equations [1] 
and [2] hold when neglecting Z)-mixing effects, which in 
the Standard Model (SM) give negligible corrections to 
7 0] and do not affect the measurement. 

Determining the angle 7 from the measurements of 
Rads and Aads requires extracting the strong phases 
by means of a Dalitz analysis of the three-body decay of 
the neutral D meson, for which the available statistics 
are insufficient. However, with the current statistics we 
can measure Rads and constrain tb hy exploiting the 
fact that in Eq. [T] |C| < 1. Since the value of is re- 
lated to the level of interference between the diagrams of 
Fig. [U high values of tb lead to a better sensitivity to 7 
in any measurement involving B D^K decays. Thus, 
Tb is a key ingredient for the extraction of 7 from other 
measurements Q. 

Both the Belle and BABAR collaborations have pub- 
lished similar measurements but in a different decay 
chain, B DK with D — Kt: 0]. Unlike those mea- 
surements, we can take advantage of the smaller value 
of ro, given by rf, = (0.214 ± 0.008 ± 0.008)% [| in 
D Kmr" decays as opposed to r|) = (0.362 ± 0.020 ± 
0.027)% [oil in D ^ Ktt decays. This implies that for a 
given error on Rads, the sensitivity to is better. 
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II. EVENT RECONSTRUCTION AND 
SELECTION 

The resuhs presented in this paper are based on 
226 xlO^ r{4:S) BB decays collected between 1999 
and 2004 with the BABm. detector at the PEP-II B 
Factory at SLAC 10]. Approximately 7% of the col- 
lected data (15.8 fb"^) have a center-of-mass (CM) en- 
ergy 40 MeV below the T(4S') resonance. These "off- 
resonance" data are used to study backgrounds from 
continuum events, e'^e~ ^ qq (q = u, <i, s, or c). The 
BaBAR detector is described elsewhere [ll|. Charged- 
particle tracking is provided by a five-layer silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH). 
In addition to providing precise position information for 
tracking, the SVT and DCH also measure the specific 
ionization (dE/dx), which is used for particle identifi- 
cation of low-momentum charged particles. At higher 
momenta (p > 0.7 GeV/c) pions and kaons are identi- 
fied by Cherenkov radiation detected in a ring-imaging 
device (DIRC). The position and energy of photons arc 
measured with an electromagnetic calorimeter (EMC) 
consisting of 6580 thallium- doped Csl crystals. These 
systems are mounted inside a 1.5T solenoidal super- 
conducting magnet. 

The event selection was developed from studies of off- 
resonance data and BB and continuum events simulated 
with Monte Carlo (MC) techniques. A large on-resonance 
data sample of B~ D^tt", — > i^~7r+7r° events was 
used to validate several aspects of the simulation and 
analysis procedure. We refer to this mode as B ^ Dir. 

Both kaon candidates are required to satisfy kaon iden- 
tification criteria, which are based on the specific ioniza- 
tion loss measured in the tracking devices and on the 
Cherenkov angles measured in the DIRC and are typi- 
cally 85% efficient, depending on momentum and polar 
angle. Misidentification rates are at the 2% level. The 
TT*^ candidates are reconstructed as pairs of photon candi- 
dates in the EMC, each with energy larger than 70 MeV 
and a lateral shower profile consistent with an electro- 
magnetic deposit. These pairs must have a total energy 
greater than 200 MeV and 118 < m^^ < 145 MeV/ c^. 
To account for the correlation between the tails in the 
distribution of the Kmr'^ invariant mass and the ir^ can- 
didate mass, we require the difference between the two 
measured masses to be within 32.5 MeV/c^ of the ex- 
pected value of TO£)o — m^o = 1729.5 MeV/c^ [l2|, retain- 
ing 90% of the signal. The remaining background from 
other B^ — > [hih2Tr^]Dhf [1| modes is reduced by re- 
moving events where the invariant mass of any ft,i/i27r'^ 
candidate, with any particle-type assignment other than 
the signal hypothesis, is consistent with the Z?" meson 
mass, retaining 92% of the signal. 

After these requirements, the background is mostly 
due to -D°- Z?" pair production in e+e" — > cc events, 
with K^TT^TT^ and D . To discriminate be- 

tween the signal and this dominant background we use 
a neural network (NNet) with six quantities that dis- 



tinguish continuum and BB events: Lq = pi and 
L2 = '^^Picos'^ 6i, both calculated in the CM frame, 
where pi is the momentum of particle i, Oi is its angle rel- 
ative to the thrust axis of the B candidate, and the sum 
runs over all tracks and clusters not used to reconstruct 
the B candidate; the angle in the CM frame between the 
thrust axes of the B candidate and of the detected re- 
mainder of the event; the polar angle of the B candidate 
in the CM frame; the distance of closest approach be- 
tween the track of the kaon candidate from the B and 
the trajectory of the reconstructed D meson (this is con- 
sistent with zero for signal events, but can be larger in 
cc events) ; the distance along the beams between the re- 
constructed vertex of the B candidate and the vertex of 
the other tracks in the event, this is consistent with zero 
for continuum events, but is sensitive to the B lifetime 
for signal events. 

The NNet is trained with simulated continuum and sig- 
nal events. We find agreement between the distributions 
of all six variables in simulation, off-resonance data, and 
B — > Dtt events. We apply a loose pre-selection on the 
NNet (0.4 < NNet < 1.0) with a 90% efficiency for sig- 
nal and a 68% rejection power for continuum, and then 
use the NNet itself in the likelihood fit described below 
to fully exploit its discriminating power. 

A i?-meson candidate is characterized by the energy- 
substituted mass mEs = \J + Po ■ PbY / Eq — p\ and 

energy difference IS.E = E^ ~ , where E and p are en- 
ergy and momentum, the asterisk denotes the CM frame, 
the subscripts and B refer to the initial e"'"e~ state and 
B candidate, respectively, and s is the square of the CM 
energy. For signal events, rriES is centered around the B 
mass with a resolution of about 2.5 MeV/c^, and AE is 
centered at zero with a resolution of 17 MeV. 

Considering both the case where the two kaons have 
the same and the opposite charge (referred to as "same- 
sign" and "opposite-sign" samples respectively), 28621 
events survive the selection described above and the loose 
requirements \AE\ < 100 MeV and mss > 5.2GeV/c2. 
While the dominant background comes from continuum 
events, there is still a non-negligible contribution from 
T(4S') BB events (denoted "BS" in the following). 
We consider separately the B Dn background, since 
it differs from the signal only in the AE distribution. For 
this decay mode the opposite-sign B~ D^it~ ampli- 
tude is suppressed by a factor « r^A^, where A « 0.22 
is the sine of the Cabibbo angle. Therefore we expect 
to find a non-negligible B — > D-k background only in the 
same-sign sample. 



III. LIKELIHOOD FIT AND RESULTS 

The signal and background yields are ex- 
tracted by maximizing the extended likeli- 
hood C = e-^'n^i A(xl)/iV!. Here N' = 
Ndk + Ncont + Nb b + No-n is the sum of the yields of the 
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FIG. 2: Likelihood fit projections of tfie NNet, /S.E, and rriES distributions separately for the same (top) and opposite (bottom) 
sign samples. To visually enhance the signal, the distributions for the latter sample are shown after cuts, with a 67% signal 
efficiency, on the ratios between the signal and the total likelihood of all the variables other than the one shown. The points 
with error bars represent the data, while the dashed, dash-dotted, and solid lines represent the contributions from continuum, 
BB, and Dn backgrounds, respectively. The dotted line represents the signal contribution, visible only in the same-sign sample. 



signal and the three background contributions (including 
both the same-sign and the opposite-sign components), 
~x = {NNet, AE, mEs}, and the likeUhood of the 
individual events (Ci) is defined as 



(*^z ) 



Nd 



K 



1 + Rads 



cont rRS 



1 + R. 



cont 



1 + Rbb 

for same-sign events and 



^^^/sKaS) + ND.fD.{^^) (5) 



NdkRaDS fWSf^-. 
NbbRbB fWS 



NcontRcont rWS ( -> \ 
-L + ncont 



1 + Rbb 



JbB i^i) 



(6) 



for opposite-sign events. In these equations we have de- 
fined R parameters for the backgrounds analogous to 
those for the signal, defined in Eq. [TJ The individual 
probability densitity functions (PDFs) / are derived from 
MC and are built as the product of one-dimensional dis- 
tributions of the three variables. The only exception is 
the TTiEs and AE PDF for the Dn background, where we 
use a two-dimensional non-parametric distribution [isj 
due to a non-negligible correlation between these two 
variables. The NNet distributions are all modeled with a 
histogram with eight bins between 0.4 and 1. The toes 
distributions are modeled with a Gaussian in the case of 
the signal, a threshold function in the case of the con- 
tinuum background, and the sum of a threshold function 
and a Gaussian function with an exponential tail in the 
case of the BB background. Finally, the AE distribu- 
tions are parametrized with the sum of two Gaussians in 
the case of the signal, an exponential in the case of the 
continuum background, and a sum of two exponentials 
in the case of the BB background. For toes and /S.E 



of the BB and continuum background, we use different 
parameters for same-sign and opposite-sign sample. 

We perform the fit by fioating the four total yields 
i^DK, Ncont, Nbb, and Ndtt), the three R variables 
and the shape parameters of the threshold function used 
to parametrize the rriES distribution for the same- and 
opposite-sign continuum background separately. Fig- 
ure [2] shows the distributions of the three variables in the 
selected sample, with the likelihood projections overlaid. 
The fit yields Rads = 0.012to°\l, Ndk = (14.7 ±0.6) x 
10^, Ncont = (239.3±2.1) X lO^, Nbb = (25.5±1.6) x 10^, 
Nd^ = (6.7 ± 0.4) X 102, R,ont = 3.05 ± 0.07, Rbb = 
0.42 ±0.07. 

Eq. [5] assumes equal efficiencies for the same- and 
opposite-sign signal samples, regardless of the difference 
in the Dalitz structure. This has been demonstrated to 
be true in MC within a relative statistical error of 4%. 
We then consider this as a systematic error on Rads- 
We also repeat the fit by varying the PDF parameters 
obtained from MC within their statistical errors and by 

tRSIWS a J i J pRSIWS 

estimatmg /co„t on ott-resonance data and Jj^jl on 
exclusively reconstructed ZJtt events. To account for the 
observed variations, we assign a 0.0076 systematic error 
on Rads- The uncertainty due to B decays with distri- 
butions similar to the signal, in particular B — > Z?(*-'7r, 
D*K, L>(*)/f*, and KK-k-k^ , is estimated by varying 
their branching fractions within their known errors and 
found to be 6 • 10~^ on RadSi and therefore negligi- 
ble. The quality of the simulation of B decays to final 
states with charm mesons that might mimic the signal 
has been checked by comparing data and MC samples in 
the sidebands of the AE distribution where these decays 
dominate. Similarly, we searched the sidebands of the 
TO^io — TOttO distribution for background from charmless 
B decays and found no evidence of it. 

Following a Bayesian approach, we extract tb by defin- 
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ing the posterior distribution 

^^^^^ ^ J p{rB,rD,£,)C{RADs{rB,rD,£,))drDd^ 

J pirB,rD,O^RADs{rB,rD,0)drDd^drB ' 

(7) 

where ^ — Ccos7, RADsifBjfDjS,) is given in Eq. [TJ 
and p{rB,rD,S.) is the prior distribution for these three 
quantities. They are considered uncorrelated, with ^ and 
rs uniformly distributed in the range of [—1, 1] and [0, 1] 
respectively. The prior distribution for rjj is a Gaussian 
consistent with r^, = (0.214 ± 0.008 ± 0.008)% [8]. The 
likelihood C{Rads) is obtained by convolving the likeli- 
hood returned by the fit with a Gaussian of width 0.0076, 
equivalent to the systematic uncertainty. 

Figure [3] shows C[Rads) and C[rB)- We set a 95% 
confidence-level (C.L.) limit by integrating the likelihood, 
starting from Rads — {tb — 0), thus excluding un- 
physical values, and we define the 68% C.L. region, for 
each variable r = Rads or ''Sj as the interval where 
C{r) > Cmtn and 68% = JcM>c^i„ ^{'r)dr. 



IV. CONCLUSIONS 



In summary, we measure the ratio of the rate 
for the B± [K^tt^t:^]dK^ decay to the fa- 

vored decay B"^ — > [K^tt'^'k^]dK^ to be Rads — 
0.012+{]:[J}g(stat)l[J;°J°(sys). While this result is consis- 
tent with and similar in sensitivity to the completely in- 
dependent previously published results |7[ , it is obtained 
using a different D decay mode. Because the measure- 
ment is not statistically significant, we set a 95% C.L. 
limit Rads < 0.039. We use this information to infer 
the ratio between the rates of the B~ — > D'^K~ and 
B- -> D^K- decays to be rs = 0.091 ± 0.059 and con- 
sequently set a limit vb < 0.19 at 95% C.L. 




FIG. 3: Likelihood function for Rads (left) and rs (right). 
The latter is obtained in a bayesian approach, assuming flat 
prior distributions for rs and ^ = CC0S7. The 68% and 95% 
regions are shown in dark and light shading respectively. 
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